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Abstract 

We present a general procedure for calculating one-loop “Casimir” energy densities for a scalar 
field coupled to a fixed potential in renormalized quantum field theory. We implement direct sub¬ 
traction of counterterms computed precisely in dimensional regularization with a definite renor¬ 
malization scheme. Our procedure allows us to test quantum field theory energy conditions in the 
presence of background potentials spherically symmetric in some dimensions and independent of 
others. We explicitly calculate the energy density for several examples. For a square barrier, we 
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divergent positive contribution near the barrier on the inside. We also carry out calculations with 
exactly solvable sech 2 potentials, which arise in the study of solitons and domain walls. 
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I. INTRODUCTION 


The weak energy condition of general relativity, the requirement that there be no negative 
energy densities, is sufficient to prevent the appearance of exotic features such as compactly 
generated closed timelike curves ilj and superluminal travel |2]. Quantum field theory 
appears to violate this condition, however. One example is the standard Casimir system of 
parallel plates, for which there is a negative energy density between the plates. However, 
this is an idealized system, which assumes a perfect conductor with an infinitely sharp and 
flat edge. A real material will have a rough surface at the atomic scale and will also appear 
transparent to very high energy modes. Since the Casimir energy is a sum over all energies, 
it will always include modes for which these effects are relevant, so this idealization could 
affect the value of the sum. In addition, since the boundary condition is imposed externally, 
there is no measure of the energy that would be required to maintain it. 1 

Although there are other ways to produce negative energy densities, for example a super¬ 
position of states with zero and two photons, these cases are constrained by averaged energy 
conditions, which require the energy to be positive when averaged along an entire geodesic. 
They are also constrained by quantum inequalities Jj, which limit the total negative en¬ 
ergy that can exist when averaging over a certain period of time. Thus it is important to 
understand problems of the Casimir type if we want to know whether quantum field theory 
protects general relativity against negative energies. 

In this paper, we reconsider the question of the energy density in such systems. To 
avoid the subtleties associated with the Casimir problem, we consider a quantum field in 
the presence of a background potential (i.e., a field with a mass that depends on position). 
In such an approach, one can choose a potential that depends only on one spatial dimension, 
and simulate the parallel plates in the Casimir problem. By generalizing the approach of 
ll, 0, {§] to local densities, we can precisely cancel the divergences in the calculation in a 
definite renormalization scheme. 

Ref. Q considered similar problems for the special case of reflectionlcss potentials, such 
as the potentials for the supersymmetric kink and sine-Gordon models. Here we present a 
general approach suitable for numerical computation, in addition to analytic calculations in 
exactly solvable models 10]. These techniques are also useful for the study of Casimir forces 


and stresses 11 


II. A SIMPLE MODEL 


To illustrate our method, we will first consider a simple model. We take a real, massless 
scalar field in 2+1 dimensions in the background of a repulsive potential V that depends on 
one spatial dimension but not the other. 

We start with the Hamiltonian density 


n 


4> 2 + (V0) 2 + U0 2 


(i) 


1 Recently, Heifer and Lang 3] showed that a frequency-independent dielectric would not be expected to 
give negative energy densities, but Sopova Q showed that negative energy densities can be achieved in 
a Casimir system with Drude-model plates, as long as the spacing is very large compared to the plasma 
wavelength. 
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and expand the field 0 in terms of small oscillations, giving 


) = 


dp 


dk 


Loc \j2n Jo 


(^k( x ) eipVe ~ Mp)ta l P + ^{xYe-^e^alJ 


( 2 ) 

where u(p) = \/fc 2 + p 2 . The J>^(x) are normal mode wave functions, which can be taken to 
be real. The sum is over the symmetric mode, +£ (x) , and the antisymmetric mode, 0^7 ( x ). 
They satisfy 

- + V(x)ip%(x) = k 2 i/i%(x) (3) 

where prime denotes differentiation with respect to x, with the normalization 


/ dx^{x)J} k ,(x) = 0 and ^ / dx^{x)^,{x) — 2nS(k - k') , (4) 

J —DO X = ~\ ^ ~°° 

which gives il)^°\x) = cos kx and ip^°\x) — sin kx as the solutions in the free case, 
V{x) = 0. 

The energy density is then 
dp 


< W ) = i-00 /o 2 mx{p) ^ \ ^ P 1 + v {x))^l{xf + ^ [xf] 

X “t“ j 


We write 




and then use Eq. © to obtain 

ro ° dk 


m - r £ 


E 


-oo 2?r Jo 27 Tu(p) x £ + _ _ 


u(p) 2 ^(x) 2 + ^-^ {j’k(x) 2 ) 


(5) 


( 6 ) 


(7) 


The integral is highly divergent, but by using dimensional regularization and introducing 
counterterms into the Hamiltonian, as discussed in Sec. nm we can render the integral finite. 
We then integrate out the transverse modes, giving 



where n is the number of transverse dimensions, later to be set to 1. The Hamiltonian 
has been renormalized by subtracting a constant term for the vacuum energy and a term 
proportional to the potential, which is sufficient to render it finite in 2+1 dimensions. 

We can relate the norms of the mode wavefunctions to the Green’s function by (see for 
example 12]) 

V'fc (x) 2 — 2klmG{x, x, k) , (9) 

x=+>- 

where G(x,x r , k ) is the Green’s function, which satisfies 

— G"(x , x 1 , k) + V(x)G(x, x', k) — k 2 G(x, x', k) = 5(x — x') (10) 
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and has only outgoing waves (~ e lfc l x l) at infinity. 

The Green’s function has the symmetry property G{x,x', k) = G{x,x', —k*)*, so for real 
k we can write 

^{x) 2 — — G(x, x, k) — 7 G{x,x,—k). (11) 

x=+- 

Tlius we can compute the energy by extending the range of integration to — oo, and using 
G^°\x, x, k ) = i/{2k), 


<W> 


ren 


p ( _ n+1 \ /*oo 

V 2 ) I , .71+1 

,, '^+3 / U 

(47t) 2 y.oo 



G(a:, x, k ) 


y(x) 

2fc 2 


+ 


n + 1 k d 2 
2 iu; 2 dx 2 


G{x, x, k) 


dk , 

( 12 ) 


where, since we are taking the massless limit, uj = \/k 2 . 

Next we would like to convert this expression into a contour integral by closing the contour 
at infinity in the upper half plane. The contour at infinity does not contribute, because for 
large, positive Im k , 

yG(i,i:,fc)^l + ^ + 0(r 4 ). (13) 

Singularities in the Green’s function in the upper half plane correspond to normalizable 
eigenfunctions of the Hamiltonian, which represent bound states. Since the Hamiltonian is 
Hermitian, the bound states must have real energies, so the singularities must lie on the 
imaginary axis and have Im k < n where /i is the mass. In this example, we have a repulsive 
potential and a massless particle, each of which is sufficient to ensure that there are no 
bound states at all. Thus the Green’s function has no singularities for Im k > 0, and the 
only contribution to the integral comes from the branch cut along the positive imaginary 
axis coming from u; n+1 . Integrating around the branch cut and using Eqs. © and © 
below gives 


{H) re n = 


2(4tt) 2 T (s±5) Jo 
and then setting n — 1 gives 


(7~i)ren — o / 

OVT .In 


K 


n+1 


„ V(x) n+1 d 2 

2 kG(x, x, k) - 1 + H-^7 —+Z+G{ X , x, k) 


2k dx 2 


2k 3 G(x, x, in) — k 2 + 


V{x) d 2 . . 

—2- k—G{x,x,ik) 


dn 

(14) 

(15) 


Once one has computed the Green’s function, this integral is straightforward, though it may 
be necessary to resort to numerical techniques. We show the calculation for some example 
potentials in Sec. IIVI below. 


III. CALCULATIONAL METHOD 
A. Model 

We will now consider the more general case of a real scalar field of mass // in the back¬ 
ground of a potential that is spherically symmetric in m nontrivial coordinates, which we 
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label by x, and independent of the remaining n trivial coordinates, which we label by y. 
The energy density is 


H = - [r + (V0) 2 + V(r)(j> 2 + p J 0 2 )= - + - V 2 (0 2 ) - 0V 2 0 + R(r)0 2 + p 2 0 2 , 


where r = |x|. Decomposing the quantum held 0 in terms of modes gives 
27rf f d n p 1 


(16) 


0(r,D,t) = 

i,iz 


T (f) J (2tt)"/2 ^ 


X 


+ 


E 




m f'O'l Jpy 


XV 


dk 


V nuj (p) 


+ ^(r)r£(D) e ™ e — ( ^a£fl7) 


where cn(p) = i/fc 2 + p 2 + /i 2 , the bound state energies are cnj(p) = y/p 2 + p 2 — /^ 2 , and the 
sum over £ gives the partial wave expansion in the m nontrivial dimensions. 

The degeneracy factor D” 1 in each partial wave is given by the dimension of the space of 
symmetric tensors with £ indices, each running from 1 to m, with all traces (contractions) 
removed 0- By the symmetry of the indices, this dimension is given by the number of 
ways to make £ indices out of 0 or more l’s, 0 or more 2’s, and so on, which is the number of 
distinct ways to place m — 1 dividers into m + £ — 1 slots. Removing all the traces requires 
subtracting the same quantity with £ replaced by £ — 2. We thus obtain 


tm (m + £— 1)! (m + £ — 3)! T(m + £ — 2) 

' £ = £\{m — 1)! ~ (£ — 2)!(m — 1)! = T(m- 1)T(£ + 1) 


(18) 


The wavefunctions ^(r) are the eigenstates of the time-independent radial Schrodinger 
equation 



m — 1 d 
r dr 


+ 


+ m — 2 ) 


+ V(r) ) i l) c t (r) = k 2 il>‘ h (r), 


(19) 


which in general comprise both bound and scattering states. The wavefunctions and creation 
and annihilation operators are normalized as follows. For the spherical harmonics, 


J Y,z(siyY r r t ,(si)dsi 


0(JL'8>i z i’ z i 


( 20 ) 


for continuum states, 


pOO 


T(f) 


r m l^( r )*^,(r)dr 




V^k p ? ^k r p' 


and for bound states 
2vrf r°° 


r(f) 


r m 1 'ip e j (r)'ip £ ,(r) dr 


>o 


n 8(k — k') 

a k ',p] = S(k~ k')8(p - p')£«'&,4 , ( 21 ) 


Sjj'Sw 
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( 22 ) 


K>'. <>■'] = [<*><4] = o K;A*;A = 

Using these expressions, we obtain the vacuum expectation value of the Hamiltonian, 

d n p 


<«> = \Y. d ? f 

n J 



E (P) I'( r ) 1 12 + f 7 W (p) I ^k (r) \ 2 

r °°dk 1 


'0 


j Vj(p) 


Wj{r)\ + 


i o 


7t o;(p) 


i^wr 


(23) 


where D 2 = + 


m— 1 d 
r dr 


is the radial Laplacian. 


B. Renormalization with one subtraction 


For positive integer m and n, this quantity diverges, as we expect since we have not yet 
included the contribution of the counterterms. Therefore we will calculate the result using 
analytic continuation in m and n from values where it is convergent. After introducing 
counterterms also depending on m and n, we will then let the dimensions go to their physical 
values while holding the renormalization conditions fixed. 

The first counterterm we will introduce renormalizes the cosmological constant. It is 
simply an overall constant in the Hamiltonian, and is fixed by the renormalization condition 
that the energy density of the trivial background V (r) = 0 is zero. The free wavefunctions 
are given by 


4>l (Cl) (r) 



r (?) 


2n 2 r 2 


-l 


Jf+t-i(kr ), 


(24) 


which, by the Bessel function identity 


\ ^ (2 q + 2i)T(2q + 1) ( _ r(2g + 1) / z\ 2c t 

4^ r(£ +1) q+e{z) — r(g +1) 2 \2/ 


satisfy the completeness relation 


^Bf|^ (0 >(r)| 2 

t 


1 k m ~ l 
(4 7 r)f- 1 2T(f)’ 


(25) 


(26) 


independent of r. Subtracting the energy in the trivial background we have 


(«} - <«} o = \ E / 

f! J 


d n p 

(2tt)" 


L 3 


E^(p)W(oi 2 + / n ^(p) (iv4(Of 2 -1^ (0) wi 2 ) 

00 dA; 1 


E 




l^(r)| 2 + 


7T £u(p) 


l«M|- 


(27) 


By completeness, in each partial wave we have 

dk 


El^)| 2 + / ^(l^(r)| 2 -|^ ( 0 ) (r)r)= 0 , 


(28) 
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so we can implement a local version of the Levinson’s theorem subtraction used in |y], giving 

d n p 


(«) - (W)o = \ Y. D ? I 


(2ir)' 


- Vp 2 + V 2 )\rf(r)\' 


dk 


+ / —(Mp) - Vp 2 + p 2 ) (iv4(r)l 2 - IV4 (0) MI 2 ) 

.In 71 " 


i« 2 E 


uAp) 


W(r) | 2 + 


dk 1 

7T u;(p) 


I^MI 1 


(29) 


This subtraction is necessary to avoid the appearance of spurious infrared singularities in 
calculations in one space dimension. These singularities also appear in dimensions less than 
one, which we will need to consider as part of the dimensional regularization process. 

Next we carry out the p integral, using 


d n p 

(27r/ 


Vp 2 + q 2 ) = 


T (f) (47r) n / 2 J 0 


P" 


l dp f a/p 2 + qA = 


£ (-^) j 

■f) 


n+a 


\n /2 


, (30) 


where we have done the integral by analytic continuation from values of a and n where it 
converges. We thus obtain 


r ( rc+n 

(W)-(W)o = - J. 


+ 


2(47t)^ 
3 dk 


- » n+1 ) W(r)| s 


7T 




+ ! / 4 D 2 (EK) n " I W( T ')l 2 + 


dk 


7T 


“ WMI' 


(31) 


where a; = £fc 2 + /i 2 for the scattering states and = \JP 2 ~ for the bound states with 
k = i/cj. 

Next, we must include the contribution of the counterterm proportional to V(r), which 
is introduced to cancel the tadpole graph. In dimensional regularization, the contribution 
to the Hamiltonian from this counterterm is 


p ( 1—n—rn \ 

nj _ V 2 / m+n— 1 

_ - , . . m+n+l P 

2(47t)— 2 — 


V(r), 


so by using 


dk T ( m ~ 2 

— (uu n+1 - n n+1 )k m ~ 3 = n™* 11 - 1 { 2 


)r(i^) 


7T 


and Eq. we have 

r(-=^i) 


H 1 = 


dk 


2(47t)^ 


- 


UJ 


n+1 n+1 


2rf(-2±l) 


V(r) 


(32) 


(33) 


M )(2_m)_ 2^ W ’‘ (r 


(34) 
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so that 


(H) ren = (H) - (H )o - (Hi) 

1 2j E D ? 


2(471)^ 




,ftn+1 ,,n+lM„//^M2 

U 


^ +1 )|^(r)| 


L 1 


+ 


<i h^" +1 - A*" +1 ) ( IV4WI 2 - IVi (0) M| 2 (l + (2 - m) V(r 


7r 


2fc 2 


+( ek)” _i i^wi 2 +r f"”- v*mi : 


We will use 

l'0fc( r )| 2 = 2&;ImG^(r, r, /c), 
where the Green’s function is defined by 


D'iG, 


’{r, r', k ) + (r) + ^ ^ £^(r, r ' 5 k) = (r — r) 


( 35 ) 


(36) 


(37) 


with the boundary conditions that it is regular at the origin and and has only outgoing 
waves (~ e lkr ) at infinity. Using G(r,r, —k ) = G(r,r, k)*, we can rewrite Eq. (ET51) as 


(H) ren = (H) - (H)o - (Hi) 


2(47t)^ 

+ r -<» 

7T 




n+1 


L 7 

n+1 


k 


H n+1 ) ‘-r ( G e (r, r, k ) - G^ 0) (r, r, fc) I 1 + (2 - m) 


y(r) 

w 


+ UJ-U 2 EG)"-GWI 2 + f°° -^-Gt^k) 

4 \ 7 T X 


(38) 


These subtractions are sufficient to render the theory finite for m + n < 3. However, it 
appears that if we set n = 1, the gamma function will cause Eq. (13%1) to diverge. In fact, as 
we will see in Appendix [0 it does not, because the quantity in brackets vanishes. But here 
we will keep n general and instead close the contour of integration in the upper half k plane. 
For sufficiently small n the contour at infinity does not contribute. There is a pole for each 
bound state at k = iEj, where Ej < ji is the bound state energy, and the contributions from 
these poles exactly cancel the sum over bound states |.16| in Eq. (1HH1) . Thus the final result 
is just given by the contribution from the branch cut along the imaginary axis from fi to oo 
resulting from a; n+1 , which contributes 

n n+ 1 (H +1 - H) n+1 ) = 2 iQ n+1 sin ^ n + 1 ^ ? (39) 

where f1 = \J k 2 — fi 2 and k = Ik. Then using the identity 

7r 


Sill 7 TZ = 


r(z + i)r (-z) 


(40) 
















we have 


(H) ren = (H) - (H )o - (Hi) 

— _1_V D m 

2(4vr)^r(^) V 


dn2nQ n+1 


Gz(r, r, in) 


Gf\r, r, in) 


x ( 1 — (2 — m) 


V(r)\ n +1 2 


2k 2 


40 2 


D;G £ (r,r,in) 


(41) 


We can now put in integer values of m and n without any divergence, as long as m + n < 3. 


Eq. (1411) can also be efficiently evaluated numerically 11 


C. Higher subtractions 


When we have m + n = 3 space dimensions, we will need to introduce a second countert¬ 
erm, | cV{x ) 2 . The first subtraction is particularly easy to define because there is a natural 
scheme, specified by the complete cancellation of the tadpole graph. Higher subtractions 
require a definition in terms of a renormalization scale, which can be chosen arbitrarily. 
In choosing this scale, we must be able to relate it to physical inputs, such as masses and 
coupling constants, in order to define a predictive theory. 

To define the counterterm precisely, we consider the two-point function n(p 2 ) in dimen¬ 
sional regularization. It diverges as we approach the physical dimension. The divergence 
is canceled by the contribution of the counterterm to the two-point function, which is just 
c. We define the renormalization scale M by taking c = —U(M 2 ). With this definition, we 
have 


c 



1 


/ \ n+m 

2(47r)^ _ T 


dE d n+m q 1 

2n (2n) n+m (E 2 — q 2 — /i 2 + M 2 A(1 — A) + ie) 2 

l*oo qn-\-m—l 

(2±=) Jo ^(4w 2 - M' 2 ) dq ’ 


(42) 


where q is the total momentum and we have integrated over the Feynman parameter A 
and the loop energy E. Typically we will choose M 2 = ji 2 , except in the case of massless 
theories, where to avoid infrared singularities we will choose a spacelike renormalization 
point M 2 < 0. 

This regulated expression is defined precisely as an analytic function of the dimension. 
Our goal is now to rewrite it in a way that allows us to incorporate it into our expression 
for the energy, Eq. m, which is also given as an analytic function of the dimension. We 
express Eq. (ESI) in terms of Green’s functions and analytically continue to express Eq. flU 
as 


c = 


ir +1 /(/yM) 


K 


m— 1 


(47r)f“ 1 2T(f) 


dhi 


2(47r)^r (2±2) J, 

/ Sr +1 /( K ,M)2 K Gf>(r,r,; K )<f K , 


(43) 


where /(/c, M) is given in terms of the hypergeometric function as 

2{m — 4)(m — 2)(2 k) 2_T71 /1 0 m3 M 2 \ 

nK ’ ’ (4« 2 -M 2 ) 3 -f sin(^f) 2 X V2’ 2’2’M 2 -4 k 2 J 


(44) 
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as we show in Appendix 1X1 

Since we will eventually take the limit where m becomes an integer, we note that 


/(«,M) = 


12 k 2 - M 2 
2k 2 (4k 2 — M 2 ) 2 


/(k, M) — 


7tk 2 (4k 2 — M 2 ) 


1 + 


4k 2 arctan ,, j) 7 , r9 

\/4k 2 -M 2 


MV4K 2 - M 2 
/(k,M) = 


2k 2 (4k 2 — M 2 ) 


for m = 1, 

for m ^ 2 , and 

for m = 3. (45) 


Eq. P^l) is now in a form where we can include it under the integral sign in Eq. m and 
obtain 


C H)ren = (H) - (H ) 0 - (Hi) - (H 2 > 

1 


2(47T)^r(^) 
V{r 


___ /*c 

E D ? / 

n J LI 


dn 2k 


fl n+1 ( G^(r, r, ?k) — G® (r, r, ?k) 


x 1- 


2k 2 


(2 — m) + E(t)7(k,M) 


n + 1, 




. (46) 


Before we can take the limit where m + n = 3, however, there is one more potential diver¬ 
gence in Eq. dUD- Our subtraction has cancelled the terms of order 1, V(r)/n 2 and V{r)/ k 4 
in the large-K expansion of the norm of the wavefunctions. But there could also be a term of 
order D 2 V (r)/ k 4 , which will generate a divergence in this case. In the renormalization of the 
composite operator T^ u , we have a renormalization counterterm 4 (^/A ~ \d x ) 4> 2 

fioj ]. Since we are considering just T 0 o here, this counterterm becomes exactly 

the form needed to cancel the remaining divergence. 2 We fix this counterterm by subtract¬ 
ing the tadpole diagram with the composite operator carrying momentum p 2 = M' 2 . Aside 
from this change, it is analogous to the tadpole subtraction above, with Vt?V(r) replaced 
by D 2 V(r). The scale M’ is then specified through the renormalization condition on the 
composite operator (and would typically be chosen equal to M ). As with M, a massless 
theory will require spacclikc M' 2 < 0, while in a massive theory we may set M’ = p. Thus 
we obtain the contribution 


(Kv) = 


(m + n — 1) 


4(n + 1) 2(47 t) 


n+1 

2 


r ( 


z / t 


dn n n - 4 2 KGf (r, r, in) D;V ^ 


M' 


■(2 —m). 


(47) 

This term is a total derivative, so it does not contribute to the total energy. We can split 
the contribution of this term beween the bulk and derivative terms so that it renders them 
both separately finite at integer dimensions, giving 


(«}«„ = 


2(4 n 


B+lp 

2 1 




/ 

n J U 


dn 2k 


k2 n+1 ( Ge(r, r, in) — G^f 1 (r, r, in) 


2 If we had chosen conformal instead of minimal coupling for the fields, which corresponds to adding the 
extra term \ ^m+n) (diid u —g IJ , v d\d > ')<j) 2 in the original Lagrangian, the divergent term would have cancelled 
automatically between the bulk term and the surface term and no renormalization would be necessary. 
However, conformally coupled theories have classical violations of the energy conditions f!3l . 
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X 


1 - XQ(2-m) + V(r) 2 f(K,M) - D r V ( r ) 


2 K 2 


8(k 2 - M rl )9? 


(2 — mf 


n + 1 


n 


” 1 1 D 2 T G t (r,r,iK) + —~ m)G?\r,r,in) 


( 48 ) 


By Eq. (1261) . the sum over i of the free Green’s function weighted by the degeneracy factor 
is independent of r, so we can pull the derivative outside in the last line, giving 


(^)ren = 


2(4vr 

xfl- 


_J_ V' D m f c 

V(r) 


dn 2 k 


2 K 2 


■(2 -m) + V{ryf{K,M) 


h2 n+1 ( G^(r, r, in) — ( r , r, in) 

D 2 V (r) 


n +1 n _! 2 


Q^D 2 G t (r,r,iK) + 


8(k 2 - M ,2 )Q 2 
V(r) , n _^ ( o). 


(2 — m)‘ 


2(k 2 - M' ) 


(2 — m)Gi’{r, r, in) 


(49) 


In this form, we see explicitly that the subtraction has cancelled the leading terms for 
large k: In the surface term, inside the derivative we have implemented the same subtraction 
of the leading behavior of the Green’s function as we found for the tadpole graph in the bulk 
term; in the bulk term, we have subtracted the leading term proportional to D 2 V(r) (which 
is derived explicitly for m — 1 in Eq. (IBlljl b For the purposes of calculation, however, it is 
often easier to work with the combined expression 


<«> 


ren 


X 


1 D m f 

(47r)^r (2±3) “ *4 
(2 - m)V(r) 




n +1 


n + 1 n 

1-— D 2 

An 2 r 


G 


e(r, r,iK) 


G 


(o) 


1 - 


, ( 2-m)D 2 V{r ) 


2k 2 


+ V (r) 2 /(ft, M) + 


8 (k 2 - M ,2 )Q 2 


(m + n — 1) 


(r, r, ik ) 
d,K . (50) 


IV. EXAMPLES WITH ONE RELEVANT DIMENSION AND ONE IRRELE¬ 
VANT DIMENSION 

A. The general case 

To illustrate our method, we would like to carry out some sample calculations in the 
case of m = 1 and n — 1. Since we are in 2 + 1 dimensions, we only need one subtraction, 
proportional to V(x). For rri — 1, the sum over partial waves reduces to a sum over the 
symmetric and antisymmetric channels. The free wavefunctions thus become 



= cos kx tpi, (°)(x) = sin kx, 

(51) 

so that 

w (0) wr + iv>h 0 ) wi 2 = i- 

(52) 

We can 

sum over the two modes to get the overall Green’s function, 



G(x, x, k ) = G + (x, x, k ) + G~(x, x, k ), 

(53) 

with 

G {0 \x,x,k ) = ^ . 

(54) 
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Then from Eq. m we have 
(V) ren = 


2(47r) 2 r ( z ^) 4 
and for n = 1 we have 


dnfl n+1 


Y( r n -|— 

2 nG(r, r, in) — 1 + - ^^-nD^G^r, r, in) 


2 n 2 20 2 


l T Ef-rf n d 2 

(‘H)ren = —— / dn fl 2 2nG(r, r, in) — 1 + —--— —-j—rG(r,r,in) 


8n 


2 n 2 fl 2 dx 2 


(55) 


(56) 


which reduces to Eq. m when n = 0. For simplicity, we will consider the massless case for 
the remainder of this section. 


B. Outside a potential with compact support 

We next consider a potential that vanishes for all \x\> a and calculate the energy density 
in this region. In this case, the only counterterm is the vacuum energy, and we get 


(^)ren = 


8n 


dn 


d 2 

n 2 (2 nG(x, x, n) — 1) — n—^G(x, x, n) 

dx 2 


The Green’s function for x, x' > a is 


G(x,x', k) = G- ( e~ ikx < + r(k)e ikx <) e* 


ikx^ 


(57) 


(58) 


where r{k) is the reflection amplitude. Thus 

2nG(x, x, in) = 1 + r{in)e~ 2hX 

and 

(H) ien = — / dnn 2 r{in)e~ 2KX . 

8tt Jo 

In the large x limit, only small n contribute in the integral. As a result, the integral depends 
only on r(0) = —1 (at k — 0 we always have perfect reflection 3 ), so we can approximate 


(59) 

(60) 


<«}. 


8n 


dn n e 


2„—2kx 


32nx 3 


(61) 


C. Square barrier 

Next we consider a square barrier with V = Vo for |x| < a and V = 0 otherwise. In 
this case, we can compute the normal mode wave functions in closed form, but must do a 
numerical integration at the end. Outside the barrier, the energy is given by Eq. (EOD with 


r = 


Voe 2na tanh 2 n'a 


2nn! + (n 2 + n' 2 ) tanh 2n'a 


(62) 


3 The only exceptions to this rule are potentials with a bound state precisely at threshold EE -LlI], which 
include reflectionless potentials. 
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and k' 2 = k 2 + Vo. 

Thus, outside the barrier, the energy is 


Vo K 2 e 2k V tanh 2k' a 

_ / (i k _ 

8n J 0 2kk! + (k 2 + k' 2 ) tanh 2k,' a 
Vq q 2 e- 2 Vy- * 1 ) tanh 2q' 

87ra J 0 ^ 2qq' + (q 2 + q' 2 ) tanh 2 q' ’ 


(63) 


where we have defined the dimensionless quantites y = x/a, q = Ka, and q' = y/ q 2 + v = n'a 
where v = Vo a 2 . Note that the integrand cannot be less than 0, so the energy outside the 
barrier is always negative. 

Far from the potential, specifically where y — 1 > l/y/v and y — 1 1/d, the contribution 

comes primarily from q -C y/v, and thus q' ~ y/v. The integral is then 


dqq 2 e 2q(y ^ 


4 (y — 1) 3 d 


(64) 


and 


<W>. 


1 

327r(x — a) 3 


(65) 


in agreement with Eq. (En¬ 
close to the potential, specifically when y — 1 -C 1/a/d and y — 1 -C 1, the contribution 
comes mostly from q y/v and g > 1. Thus q' ~ q and tanh 2q' ks 1. The integral becomes 


1 

8(2/ -1) 


and 


Inside the barrier, we have 


<«>. 


Vo 

647r(a; — a) 


( 66 ) 

(67) 


, i {k' cos k'(x < + a) — ik sin k'{x < + a))(k' cos k'(x > — a) + ik sin k\x > 

X,X ’ k' 2kk' cos 2k'a — i{k 2 + k' 2 ) sin 2k'a 

where k! = y/k 2 — Vo, so we can write 

1 (k 2 + hi' 2 ) cosh 2 n'a + 2kk' sinh 2k' a + Vo cosh 2k' x 
X ’ Xl ' 2 k! 2kk! cosh 2k! a + ( k 2 + k' 2 ) sinh 2k' a 

We can then split the energy into two parts, 


a)) 

( 68 ) 

(69) 


{'H) ren — Eq + E\(x) , (70) 

where £7 depends on position, but E 0 does not. 

The position-independent part is 

E f°° ^ f r3 ( r2 + r/2 ) cos ^ + 2 ac/c' s i n h 2/da 2 Vo | 

0 87t 7 0 i k! 2kk' cosh 2k' a + (k 2 + k' 2 ) sinh 2k' a 2 / 
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V n 2 f , 1 2k + k tanh 2k a 

87 r J o 2k 2 kk' + ( k 2 + k' 2 ) tanh 2k' a 
_Vla r°° d 1 2q + q' tanh 2q' 

8n J 0 ^ 2q' 2qq' + (q 2 + q' 2 ) tanh 2 g' 

and is always negative. In the limit where v 1, we can approximate tanh 2q' fa 1 to get 


E 0 = 


y3/2 
[0 _ 

12?r ‘ 


(72) 


The position-dependent part is 


EAx) = 


1 

87 r 


d,K 


K 


(! 2k' 2 — k 2 )VqCos\\2k'x 


Vr 


0 


87ra 


k! 2kk' cosh 2k' a + (k 2 + k' 2 ) sinli 2k' a 
q dq (2 q' 2 — q 2 ) cosh 2 q'y 

q' 2qq' cosh 2 q' + ( q 2 + q' 2 ) sinh 2 q' 


(73) 


and is always positive. 

Note that the the dominant term in the integrand in Eq. CUD is suppressed by e 2?, d y ) < 
e — 2 \A7(i— 2 /) xhus, far from the edge of the potential, where 1 — y ^-/y/v, E\ is negligible. 

Close to the edge of the potential, with 1 — y -C 1/y/v and 1 — y 1 , the integral is 
1 /( 8(1 — 2 /)) and 

EAx) w- -4- -r (74) 

647r(a — x) 

which cancels, in a principal value sense, the divergence outside the barrier. 

The sign of the energy density at the center of the barrier depends on the competition 
between the position-dependent and position-independent parts. For large v, the position- 
dependent part is suppressed in the center, and the energy density is negative. For small v, 
it is positive. The total energy density is shown for several values of v in Fig. |H 

It has long been known that the energy density near a perfectly reflecting boundary is 
zero if one uses the “conformally coupled” stress-energy tensor, but diverges if one uses the 

argue 


minimally coupled one, as we have done above. Kennedy, Critchley, and Dowker 3 


that since the total energy is the same in the two cases, there must be a surface energy 
associated with the perfect conductor in the minimal case. Ford and Svaiter 0 found that 
the surface energy could be seen by allowing the boundary to fluctuate. 

Here, we can see the situation by approximating a perfect conductor by a square barrier 
with a fixed and Vo —> 00 . We can produce the conformal Hamiltonian density by including 
half the value of the total derivative term, 


77conformal = ^ V 2 (0 2 ) - ^0V 2 0 + V (r)0 2 . (75) 

Z o Z 

This choice gives zero energy outside the barrier and removes the divergence of the energy 
density everywhere inside. With the minimally coupled Hamiltonian, the energy outside 
goes to Eq. (ESI) as Vo becomes large, while the positive energy inside clusters ever closer to 
the boundary, as shown in Fig. El Since the change to the total derivative term does not 
affect the total energy, we can see that the “surface energy” located just inside the boundary 
cancels the divergent negative energy outside. 
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o /o 

FIG. 1: Energy density in units of V 0 ' for the square barrier of width 1 and heights 0.05 (dashed), 
1 (solid), and 5 (dotted). As Vo increases, the positive energy becomes concentrated more and 
more near the edge of the barrier. In units of V Q , the outside energy decreases with Vo, but in 
absolute terms it approaches a fixed limit given by Eq. ©■ 


V. EXAMPLES WITH ONE RELEVANT DIMENSION AND TWO IRRELE¬ 
VANT DIMENSIONS 


A. The general case 


To carry out calculations in 3 + 1 dimensions, we now need subtractions proportional to 
V(x), V(x) 2 , and V"(x). We will use the renormalization scheme defined in Section HTTI For 
m = 1, using Eq. (gSD and evaluating Eq. m with m = 1 and n — 2 gives 


(H) ren = 


127T 2 


dutt 3 


V(x) 2 {12k 2 - M 2 ) 

2k 2 (4k 2 - M 2 ) 2 ~~ 4(k 2 - M' 2 )Q 2 


2kG(x, x , in) 


3k d 2 

2Q 2 dx 2 


G(x, x, ik) — 1 + 


V{x) 
2 K 2 


V"(x) 


(76) 


where the Green’s function has again been summed over the symmetric and antisymmetric 
channels. Again, we will restrict our attention to massless fields for simplicity. 


B. Outside a potential with compact support 


The wave functions and Green’s functions are just as in Sec. IIVI Again, since the potential 
vanishes, the only counterterm is the vacuum energy. Thus Eq. dZHD reduces to 


(H) ren = 


1 


12tt 2 


dn 


3 0 d 2 

k 3 (2 kG(x, x , k) — 1)- k 2 -—G(x, x, k) 


2 dx 2 


outside the potential, and so from Eq. 

1 f 


(H) ren = 


67T 2 


dn K 3 r(in)e 


—2kx 


(77) 


(78) 
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In the large x limit, we can again take r{in) ~ r(0) = —1, to get 


(^)ren = 


67T 2 


= 


167r 2 a; 4 ’ 


( 79 ) 


a well-known result. 


C. Square barrier 

Outside a square barrier with width a and height Vo, the reflection coefficient is given by 
Eq. (jH2J), and the energy is 


V 0 f°° , K 3 e 2k ^ x a Hanh2 n'a 

[H) ren = / «K 


6n 2 Jo 2kk' + (k 2 + k' 2 ) tanh 2 k' a 

Vo f°° q 3 e ~ 2q ( y ~V tanh2g' 

dq - 

^ Onn> 


67T 2 a 2 Jo 2 qq' + (q 2 + q' 2 ) tanh 2q' ' 

Far from the potential, we approximate q' ~ y 7 v q. The integral is 

dqq 3 e- 2q ( y -V 3 


(80) 


8(y - l) 4 u 


(81) 


and 


<W>r 


167r 2 (a: — a) 4 ’ 

in agreement with Eq. tED- 

Close to the potential, we approximate q' ~ q and tanh 2q' cs 1. The integral becomes 


(82) 


1 


and 


(W)i 


16(2/ - l) 2 
C 0 


967t 2 (x — a) 2 

Inside the potential, we need the renormalized form, 


(83) 


(84) 


(H) re „ = 


1 


12vr 2 


dn 


2kG - 1 + 


V{x) _ (12 k 2 - M 2 )V(x) 2 \ _ 3 2 d 2 ' 
2k 2 2k 2 (4k 2 — M 2 ) 2 y 2 " da; 2 


where M 2 < 0 is the spacelike renormalization point. 

For the square barrier, we get a position-independent part, 


(85) 


E n = 


12vr 2 


1 f°° f ft 4 (ft 2 + ft /2 ) cosh 2 k' a + 2 kk' sinh 2k' a 

dn 


k! 2kk' cosh 2k' a + (k 2 + ft' 2 ) sinh 2k' a 
kVo k(12k 2 + M 2 )VJ 


-k° + 


t/2 roc 
V ° 1 dK 


127T 2 


'0 


2 2 (4k 2 + M 2 ) 2 

k 2k + k! tanh 2k' a k(12k 2 + M 2 ) 


2k' 2kk' + (k 2 + k , j ) tanh2 k 'a 2(4k 2 + M 2 ) 2 
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Vq f°° f q 2q + q' tanh 2q' q(l2q 2 + t 2 )'\ 

12n 2 J 0 q l 2q' 2qq' + (q 2 + q' 2 ) tanh 2 g' 2(4g 2 + t 2 ) 2 / ’ 

where M 2 = —M 2 and t = Ma. 

We can isolate the dependence on the renormalization scale by using 

f,f 3 g(12g 2 + t 2 ) } 3 , t 2 1 

X M S(/ 2(4<2 2 + f 2 ) 2 / 16 « + 8 


( 86 ) 


(87) 


to obtain 

p + I , f./g 2 g + ^tanhV ,3^ 

0 127T 2 yi 6 Vo 8 J 0 12 g' 2qq' + (g 2 + g /2 ) tanh 2 g' 8 g'Jy 


In the limit where 3> 1, we can approximate tanh 2q' rs 1, the integral gives —7/32, and 
we obtain 

W 1 " 


Er\ = 


V 2 

v 0 

64:/r 2 


In + 
M 2 2 


(89) 


consistent with the result obtained from the effective potential [ 2 | 
The position-dependent part is 


1 r°° _ k 2 (3k' 2 — k 2 )V 0 cosh 2k'x 

127 r 2 J 0 ' k! 2kk' cosh 2k' a + (k 2 + k' 2 ) sinli 2k' a 

Vo r q 2 (3q' 2 — q 2 ) cosh 2q'y 
127 r 2 a 2 J 0 ^ q' 2qq' cosh 2 q' + (g 2 + q' 2 ) sinli 2 q' 


(90) 


and is always positive. 

Far from the edge of the potential, E\ is negligible. Close to the edge, where we can 
approximate q' ~ q and sinli 2q r; cosh 2 q ~ e 2 g / 2 , the integral becomes 1 /( 8(1 — y ) 2 ), and 


£^i(x) R 


Vo 

967r 2 (a — x) 2 ’ 


(91) 


which cancels the divergence outside the barrier. 

These results do not reflect any contribution from the V"(x) counterterm. In this case 
it vanishes for all |x| ^ a, since the potential is constant. Furthermore, the contribution 
to the total energy from this term is also zero, since it is a total derivative. If we imagine 
that the square barrier represents the limit in which a smooth potential gets steeper and 
steeper, we will find large equal and opposite contributions to the energy localized in the 
tiny region on both sides of the boundary. As long as we average over larger distance scales, 
this contribution will always cancel out, so it can be ignored in the square barrier limit. 


D. The sech 2 potential in 3 + 1 dimensions 

Finally, we consider the potential analyzed in 2 + 1 dimensions in na, 

V(x) = c 2 sec\i 2 (x/a ), (92) 
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FIG. 2: Quantum energy density due to the potential of Eq. & for c = 1/a, n = 0, and 
M 2 = M' 2 = —1/a 2 in units where a = 1. 


which arises frequently in soliton models. It is exactly solvable in terms of associated Leg¬ 
endre functions. For c 2 a 2 = — £(£ + 1) with integer i it becomes reflectionless. The Green’s 
function at coincident points is 

G(x,x,in) = — r(l + na + s)r(Ka — s)P7 Ka (tanh(o;/a))P“ Ka (— tanh(a;/a)). (93) 

where is the associated Legendre function as defined in ^1] for — 1 < x < 1, and 

s = (\/l — 4c 2 a 2 — l)/2. Plugging this into Eq. (I7H1) . we have 


(W)ren = 


127T 2 


dn Q 3 —r(l + na + s)r(«:a — s) 
3 k d 2 


X V 2K 29? dx 2 

c 4 sech 4 (a;/a)(12K 2 — M 2 ) 
2k 2 (Ak 2 -M 2 ) 2 


(tanh(x/a))P s Ka (— tanh(x/a)) — 1 + 
c 2 sech 2 (a:/a)(3 sech 2 (x/a) — 2) 


c 2 sech 2 (a;/a) 


2k 2 


+ 


2a 2 (k 2 - M' 2 )Q 2 


(94) 


which can then be computed numerically. Figure Ogives this energy density as a function 
of x for particular values of the parameters. 


VI. CONCLUSIONS 

We have seen how to address the question of generation of negative energies through 
quantum fluctuations in the robust language of quantum held theory, where ambiguities 
associated with idealized boundary conditions are absent. This approach implements stan¬ 
dard renormalization procedures and is applicable to generic background potentials that 
are spherically symmetric in some dimensions and independent of the rest. Such potentials 
typically arise, for example, from topological defects or other extended objects. By using di¬ 
mensional regularization, we have implemented a precise renormalization scheme, using only 
local subtractions for both the first- and second-order diagrams. We expect that this general 
formalism, together with fermion scattering theory in fractional dimensions developed in Q, 
will allow these results to be extended to fermions and gauge fields. 
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In the case of the square barrier, we have recovered the negative energy associated with 
perfect reflection at large distances from the barrier, and we have seen that the divergent 
negative energy outside the barrier is canceled by positive energy immediately inside. In a 
realistic example in which one includes the energy associated with the background potential, 
such cancellations might lead the averaged null energy condition to be obeyed even though 
the weak energy condition is violated h|. Finally, we have calculated the energy density 


for a smooth background representing a domain wall in 3 + 1 dimensions. 


VII. NOTE ADDED IN PROOF 


Ref. |22j has calculated the surface tension for a bosonic 0 4 kink domain wall (and also its 
supersymmetric generalization) using an on-shell renormalization scheme, in space dimension 
one through four. In the language of the present paper, this calculation corresponds to the 
case of a = 2//i and c 2 = — 3 /i 2 /2 in the potential of Section IV Dl renormalized with 
M — //, setting m = 1 and n to zero through three. The surface tension is obtained by then 
integrating this result over the one nontrivial dimension. (The choice of M' does not affect 
this calculation because the total derivative term integrates to zero.) Using the formulae 
in the present paper to carr y o ut this calculation, we obtain results in agreement with the 
bosonic calculations in Ref. M. 
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APPENDIX A: CALCULATION OF / 

By comparing the last line of Eq. (TI21) with the first line of Eq. © , we require 


'■oo qn+m— 1 


rq?) 


on+1 ..„_i J,„ _ 1 

—. U f\j / \ L lj / . \ I / . 

T (f) r (^) Jo u ( 4u - M ) 


dq, (Al) 


where u = \J q 1 + fi 2 and D = y//t 2 — /i 2 . 

Let us change variables on the left from k to L = D 2 = n 2 — /jl 2 and on the right from q 
to L = q 2 to get 


rpJ 2 ) 


f fa+ivy„- 2 A/(*,M) ■ = 1 r 

r(f) u(4 


dL 


(A2) 


with n = sjL + /i 2 on the left and u> = \JL + /i 2 on the right. We can write 


F(1 + a) Jo 


dL L a h(L + /i 2 ) = 


r(i + p) Jo 


dL L^j (L + jj 2 


(A3) 
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with 


a = 

(n + l)/2 

(A4) 

P = 

(n + m — 2)/2 

(AS) 

h(x) = 

f(y/x, M) 
r (¥) 

(A6) 

j{x) = 

1 

(AT) 

y/x(4:X — M 2 ) ’ 


Denote the difference in the exponents as 5 = a — (3 = (3 — m)/2. If 5 is a positive integer, 
the desired relationship is just integration by parts, and 

= (-£)* (A8) 


To extend this formula to non-integer 5, we write j in terms of the hypergeometric function, 


j(x) 


—— ( --- 

20F V4 x- M 2 


d 

r (d) 2 f, 


(1 3 M 2 \ 

\2’ 2’ M 2 — Ax) ’ 


(A9) 


with d = 3/2. The operator (—d/dx) just increments d in Eq. (1 ADD . so we conjecture that 
the same relationship holds for all 5 , and thus that the desired h is given by Eq. (IA9J) with 
d = 3/2 + 5 = 3 — m/2. One can check that the conjecture is correct by explicitly performing 
the integrals in Eq. (1A3I) . which both give 


1 


2v/? 

Finally we End 


M- 


2 \ 0 - 1/2 




- P ) 2 Tj 


1 1 


M 2 


2’2 2’ M 2 — 4/i 2 


(A10) 


M+ r hi r ^-T 


m 


m \ 


(4k 2 — M 2 ) 3_m / 2 sin 


2 ) V 4k 2 - M 2 

1 

2 


3—m/2 


K 2-m 
K 2 


p , 1 o _ 3. M 2 

1 1 2’ 2 ’ 2’ M 2 -4k 2 


2(m — A)(m — 2) (2k) 2 m /I m3 
-2Ai -, 3 - —; -, 


ixm 

2 


M 2 


2 ’ 2’ M 2 -4k 2 


(All) 


APPENDIX B: LOCAL SUM RULES 


General case 


We used the analytic properties of the Green’s function as a mathematical tool, enabling 
us to carry out calculations efficiently on the imaginary axis. In so doing, we avoided the 
apparent singularity in the gamma function coefficient of Eq. USD for odd n. Nonetheless, 
this expression should be a valid result, finite for m + n < 3. As in the case of the total 
energy |d|, the quantity in brackets must vanish for n — 1 in each partial wave individually. 
Furthermore, the combination of the first two terms in brackets vanishes separately from 
the total derivative term. These cancellations depend on a local analog of the sum rules for 
the phase shift given in 23[ 24[, which we demonstrate below. 
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When similar apparent divergences arise in the calculation of the total energy, they are 
canceled according to generalizations of Levinson’s theorem [a! 0. For a system with 
spherical symmetry, in each partial wave i these sum rules take the form 


^2(~ K £j) N + f 

Jo 


k2N i I - 


N 

J2 5< t\ k ) I dk = 0, 

3=1 


(Bl) 


where the bound states have kgj = iKgj, Sg(k) is the scattering phase shift, and S^ s \k) is the 
scattering phase shift computed at order s in the Born approximation. 4 The N — 0 case gives 
Levinson’s theorem. Like Levinson’s theorem, these identities apply to general potentials in 
scattering theory and hold in each partial wave i individually. Also like Levinson’s theorem, 
they are modified for the case of the symmetric channel in one dimension, as discussed in 

0. 

We have a relationship llj, 12j between the phase shift, the change in the density of 
states, and the norm of the wavefunction, 


1 dSi(k) 
Ti dk 


Pti k ) ~ pf\ k ) 


1 

GfU 


drr 


m—1 



(B2) 


where the zero superscript indicates a quantity evaluated in the free case. This equation 
also holds order by order in the Born approximation. Using these relations we can rewrite 
Eq. m as 


2vr- 

rjf) 


N 


drrm 1 ( X}(-«<j)*l^(r)| 2 + ^ / k2N ( IV4(0| 2 - l^ (S) ( r )! 2 ) dk ) = °> 


s=0 


(B3) 

where i/jf} s \r ) is the Born approximation to the wavefunction computed at order s (the free 
wavefunction is the order zero term). The identities we need for the present application are 
simply the slightly stronger condition that Eq. (El holds for each r individually, rather 
than just as an integral. We can exploit the connection to the Green’s function that was 
used in |23,0| to prove this result as well. 

The case of = 0 is particularly simple, because we know that 


iV’jMI 2 + - [ - \ipi {0) (r)| 2 ) dk = 0 (B4) 

A J 0 


by completeness; it is just the difference between the expectation value of a constant com¬ 
puted in the free and interacting bases. (After summing over the spectrum, each term is 
independent of r.) For higher N, we would like to show that 

IV4 W M| 2 | dk (B5) 


1 /»oo / N 

£(-'4)' v ldMl 2 + - / ivtMl 2 - £ 

I n Jo \ s =0 


4 In general, these identities continue to hold even if one subtracts N' orders in the Born approximation 
for any N 1 > N. However, there are some restrictions on such oversubtractions in the symmetric channel 
in one dimension 
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is zero. 5 We employ the relationship in Eq. m between the norm of the wavefnnction and 
the Green’s function to rewrite this expression as 

G ( °\xiX,k)^j dk , (B7) 

where we have extended the integral to the entire k axis by the symmetry of the integrand. 
To show this expression is zero, we would like to do the k integral as a contour, closed 
in the upper half plane. The singularities in the full Green’s function correspond to bound 
states, and will exactly cancel the explicit contribution from the bound states 0- The Born 
approximation has no singularities (since it does not see the bound states). Thus we are left 
with the contour at infinity. However, it does not contribute because we have subtracted 
enough Born approximations to ensure that the integrand falls like 1/1A:| 2 at large \k\ j25j |. 




2 \N\j,£ 
tj/ 


mr)\ 


1 

+ - 

7r 


/•OO / N 

/ k 2N+1 Im I Gi(x , x , k) — ^ 

V s=0 


2. The symmetric channel 


In one dimension, we have to consider the symmetric channel, which can have additional 
singularities at k = 0. Such singularities, for example, lead to an extra 1/2 in Levinson’s 
theorem 0, relating the phase shift at k = 0 to the number of bound states. We have 


5s(0) = vr 



as opposed to the usual 


5(0) = 7T n . 


Analogously in our problem, Eq. (IB4I) must be modified to 


E^)i 2 


i i 
2L + n 




dk = 0, 


(B8) 

(B9) 


(BIO) 


where L is the size of the system. Subtracting A- reflects the contribution from the state 
i]){x) = const in the free spectrum. This state is “half-bound”: While any potential will have 
a k = 0 state in the symmetric channel, in this case the wavefnnction goes to a constant at 
infinity. Such states contribute to the spectrum with half the usual residue for a bound state, 
as the name indicates. (Generically a state with k = 0 will approach a line with nonzero 
slope, in which case no special treatment is necessary.) If a potential has a half-bound state, 
making the potential arbitrarily more attractive introduces a new bound state in the theory, 


5 As shown in [2j, for reflectionless potentials in one dimension there is a stronger version of the first local 
sum rule, 

I VAGI 2 J + h 2 + \Mr)\ 2 - IV>fc (0) (r)| 2 =0, (B6) 

3 J 

which reduces to Eq. iBdll when integrated over k. It might be possible to find analogous results for the 
higher sum rules as well. 
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and making it arbitrarily more repulsive eliminates the half-bound state. 6 There will be an 
analogous contribution to the energy density, so this term will cancel when we pass from 
Eq. (1771) to Eq. (171 and the rest of the derivation of the energy density is unchanged. 

For the other sum rules needed in our problem, however, we always multiply by enough 
powers of k to cancel any anomalous effects coming from states at k — 0. We would have to be 
more careful if we do additional Born “oversubtractions,” in which case we could encounter 
additional terms analogous to those found in jUJ. We can always avoid these problems as 
long as each ultraviolet Born subtraction is preceded by a corresponding infrared Levinson 
subtraction. For the first Born subtraction, the corresponding Levinson subtraction was 
done using Eq. (IB4D . Higher Levinson subtractions would use local analog of the higher sum 
rules in 61. 


3. Local subtraction 

For Casimir calculations it will be convenient to slightly modify the N — 1 sum rule. Our 
renormalization procedure subtracts not the full first Born approximation, but rather just 
a local part of it. However, this replacement does not affect the sum rule. For example, to 
apply the results of Section Qn Bl for m — n— 1, we write 

\ M x )\ 2 

x=+- 


and subtract only the term directly proportional to V(x ), rather than all terms that are first 
order in the strength of the potential. However, the additional terms, proportional to the 
derivatives of H(x), do not introduce any singularities in the integral and do not affect the 
contour at infinity because they fall like 1/k 4 or faster. Therefore, this modification does 
not affect the proof of the sum rule. This result allows us to apply the sum rule to Eq. (ESD- 
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2 k 2 
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sin 2k (y — x) + 

V"\y) 


2k 2 


8 k 4 


x ! V{x) V"{x) + /■« 


2k 2 


8 k 4 


dy 

dy 


8 k 4 

v""{y) 

16A; 5 


cos 2 k(y — x) + 
sin 2k(y — x) + 


(BH) 


6 A reflectionless potential will always have a half-bound state, because it must have Ss(k) = for all 

k. If this equality is to hold at k = 0, to reconcile Eqs. IB8I) and (IB9I) there must be a half-bound state, 
which contributes only a half to the number of bound states. The half-bound state in the free case (which 
is reflectionless) is just a consequence of this requirement. 
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4. One irrelevant dimension 


With the sum rules in hand, we can now extract a finite result from Eq. (1351) . Near n — 1 
we have 

r f— n « —-— and a n ~ l fa (1 + ——-logo 2 ) (B12) 

\ 2 / 1 2 

so that in the n —> 1 limit we obtain 


(ft) ren - 8n ^2 D . 


^(a;j) 2 l °g(u;j) 2 |^(r )| 5 


e 

+ J o ^w 2 loga; 2 ^fc(r)| 2 - «W(r)| 2 ^l + (2-m) y ^^ 

+ T^r fctog(u# 2 |$(r )| 2 + J o ^f lo & u Vk( r )\ 


2k 2 


(B13) 


The local sum rule ensures that the scale of the logarithm does not affect the final result. 
In addition, the limit p ^ 0 is smooth (except when n = 0 and m — 1, where we have 
the usual infrared divergences of one-dimensional field theory). If we extend the range of 
integration in Ea. (IB 131) as in Ea. m, and then close the contour in the upper half plane, 
the branch cut associated with logcu 2 will reproduce Eq. dm 
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